Observation of magnetocoriolis waves in a liquid metal Taylor-Couette experiment 
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The first observation of fast and slow magnetocoriolis (MC) waves in a laboratory experiment is reported. 
Rotating nonaxisymmetric modes arising from a magnetized turbulent Taylor-Couette flow of liquid metal are 
identified as the fast and slow MC waves by the dependence of the rotation frequency on the applied field 
strength. The observed slow MC wave is damped but the observation provides a means for predicting the onset 
of the Magnetorotational Instability. 
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Hydromagnetic waves are a ubiquitous feature of both geo- 
physical and solar dynamo models as well as models of astro- 
physical accretion disks. Such waves were first observed in 
experiments [1] with liquid metals using a sufficiently strong 
magnetic field that the Lorentz force could act like tension 
on a string and support Alfien waves [2]. When the liquid 
is in rapid rotation, these waves are modified by the Corio- 
lis force. The resultant magnetocoriolis waves [3] are a hy- 
brid of Alfven waves and inertial waves [4] . Magnetocoriolis 
waves (MC waves) are used to explain the secular variation 
of the Earth's magnetic field over the course of hundreds of 
years [5] and the redistribution of angular momentum in the 
Sun [6]. They are a special case of the more general Mag- 
netic Archimedes Coriolis (MAC) waves from dynamo the- 
ory [7, 8]. 

Despite the importance of MC waves in rotating conduct- 
ing fluids and plasmas there is scant experimental evidence 
of their existence and of their relationship to various im- 
portant astrophysical phenomena such as the dynamo or the 
magnetorotational instability (MRI) [9-11]. Recent experi- 
ments [12] on a liquid sodium spherical Couette device and 
simulations [13] have found several different types of hydro- 
magnetic waves, but there is ambiguity about their identifica- 
tion as MC waves. A similar experiment in Maryland found 
inertial waves [14], but the applied fields were too weak to 
observe Lorentz force effects. 

In this Letter, we report the first clear identification of the 
combined fast and slow MC waves in a laboratory experi- 
ment. Through measurements of the radial magnetic field in 
a liquid metal Taylor-Couette flow, we observe two rotating 
modes that follow the rotation speeds expected for the fast 
and slow MC wave. We also demonstrate through a local sta- 
bility analysis that with the addition of sufficient flow shear, 
the slow MC wave can be destabilized to produce the MRI. 
Using the observed frequencies of the waves we infer from 
the local dispersion relation that the modes are damped and 
obtain a method of determining the threshold for the MRI. 

The Princeton MRI experiment is designed to study the sta- 
bility of a rotating sheared flow of liquid metal with an ap- 



Experimental Parameters 


Symbol 


Units 


Value 


Height 


h 


cm 


27.9 


Inner Cylinder radius 


n 


cm 


7.06 


Outer Cylinder radius 




cm 


20.30 


Density 


P 


g/cm 3 


6.36 


Kinematic viscosity 


V 


cm 2 /s 


2.98 x 10~ 3 


Magnetic diffusivity 




cm 2 /s 


2.57 x 10 3 



TABLE I: Physical parameters of the apparatus [15] and liquid 
metal [17]. 



plied magnetic field coaxial with the rotation axis. The ap- 
paratus has been described elsewhere [15] and has already 
demonstrated the ability to generate high Reynolds number 
shear flow in water with angular momentum flux comparable 
to viscous transport (a null result in trying to demonstrate sub- 
critical hydrodynamic instability [16]). The volume between 
the concentric rotating cylinders is filled with GalnSn, a gal- 
lium eutectic alloy [17]. The dimensions of the experiment 
and properties of the liquid metal are given in Tab. I. 

The equations describing the evolution of a rotating shear 
flow with a background magnetic field are given by the mag- 
netohydrodynamic (MHD) equations in a rotating frame: 
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where v and B are solenoidal fields representing the veloc- 
ity and magnetic field, CI is the angular velocity, and p is the 
pressure. The generalized pressure P incorporates the mag- 
netic and centrifugal pressure terms. 

The background field and angular velocity are given by 
Bo = Bqz and CI = £lz. If we assume harmonic perturbations 
of the velocity and magnetic field and linearize the resulting 
equations, we obtain the dispersion relation in cylindrical co- 
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ordinates [18]: 
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where Icq = m/r for integer m, cb = CO m£l is the Doppler- 
shifted complex frequency, = r]k 2 and y v = vk 2 are the 
resistive and viscous damping rates, COa = k z Bo / ^/poP is me 
Alfven frequency, and COr = (C, —2)£lk r k§ /k 2 is the Rossby 
wave frequency [19]. We have quantified the flow shear by 
introducing the vorticity parameter 



CM: 
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Note that £ is only constant if the rotation profile follows 
a power-law dependence £2(r) oc r^ -2 . For uniform rotation 
with no shear £ = 2. The Rayleigh criterion [20], which gov- 
erns axisymmetric hydrodynamic centrifugal stability of ro- 
tating shear flow, is given by £ > . 

We can gain insight into the basic waves for this dispersion 
relation by examining limiting cases. In the absence of rota- 
tion, the dispersion relation reduces to 



(co - iy v ) (co - iy^) - col = 



(3) 



which describes the damped shear Alfven wave. These waves 
have a transverse polarization due to the incompressibility of 
the fluid (in a compressible fluid there is also a longitudinal 
magnetosonic wave). The complex frequency is CO = i(y v + 



y 7] )/2±^jcD 2 -(y v -y 7] ) which shows that the Alfven wave 
is viscously and resistively damped and has a real frequency 
of ±C0a in the absence of dissipation. Variations in the flow 
along the direction of the magnetic field tend to be eliminated 
in a highly resistive fluid (such as a liquid metal) due to this 
damping. 

Assuming rotation without shear and without an applied 
magnetic field, the dispersion relation reduces to 



(co-iy v ) 2 + (2£lk z /k) 2 = 



(4) 



which describes inertial waves. Inertial waves also have a 
transverse polarization but are peculiar in that the restoring 
force, provided by the Coriolis effect, acts orthogonally to 
the displacement. The resulting motion of a displaced fluid 
element is circular precession. The complex frequency is 
CO = iy v ± 2£lk z /k so the wave is viscously damped with a real 
frequency between ±2£2. Note that there is no dependence 
on the wavelength. Akin to Alfven waves, inertial waves ho- 
mogenize the flow along the axis of rotation due to viscous 
damping, consistent with the Taylor-Proudman Theorem [4]. 

Together these waves make up the hybrid magnetocoriolis 
wave [3]. Since there are two restoring forces acting on dis- 
placed fluid elements, there are two possible situations. The 
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FIG. 1 : The real part of the dispersion relation for (a) damped Alfven 
waves as given in Eq. 3, (b) inertial waves as given in Eq 4, and (c) 
the fast (dashed red line) and slow MC waves (solid blue line), (d) 
The growth rate of the fast and slow MC waves. Due to resistive 
dissipation, the MC waves collapse to the inertial wave dispersion at 
high k. When there is sufficient flow shear, the slow MC wave be- 
comes unstable at low k. The parameters for the dispersion relation 
are Bq = 4 kG, Q. = 42 rad/s, and £ — 0.25 (maximum design param- 
eters of the apparatus) and the fluid properties are given in Tab. I. 



Lorentz and Coriolis forces may act together, stiffening the 
system and producing the higher frequency fast wave, or the 
two forces may oppose one another to produce the lower fre- 
quency slow wave. The requirement for observing strong ro- 
tational effects on the Alfven wave is rQ^fiop /Bq ~ 1. The 
Alfven frequency experiences a splitting due to the breaking 
of the degeneracy of the two roots by the presence of rota- 
tion. When the flow has sufficient shear, the slow MC wave 
becomes stationary and unstable at low k as seen in Fig. 1. 
This instability is the MRI which in addition to causing tur- 
bulent transport of angular momentum in accretion disks may 
also be related to geomagnetic jerks [21]. 

The results reported here were obtained when an axial mag- 
netic field was applied to turbulent rotating shear flow. The 
outer cylinder was kept stationary while the inner cylinder ro- 
tated. The rings at the end caps were coupled to either cylin- 
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FIG. 2: A contour plot of the radial magnetic field near the surface of 
the outer cylinder constructed from a least squares fit of data from the 
Mirnov coil array. The locations of the coils are depicted by dots and 
the contour levels are given in Gauss. The applied field was 4.3 kG. 



der in what is referred to as the "split" configuration [15]. The 
flow was hydrodynamically unstable since £ < 0. The inner 
cylinder and inner ring were set in rotation at 6.7 Hz and an 
axial magnetic field between 1.7 to 4.3 kG was applied to the 
turbulent flow. The induced radial magnetic field fluctuations 
were measured by an array of Mirnov coils[22] positioned just 
outside the outer cylinder. A nonaxisymmetric mode was ap- 
parent from the coil array as seen in the snapshot of the field 
shown in Fig. 2. The image is created by fitting the signals 
from the array to a spatial Fourier mode model using standard 
least squares fitting. From the model we obtain the Fourier 
amplitude and phase as shown in Fig. 3. The variation of the 
mode amplitudes with magnetic field strength is not linear as 
would be expected for advection by differential rotation (the 
a effect). 

The mode rotation rates are measured by calculating the lin- 
ear slope of the Fourier phase with time for each field strength. 
The results are shown in Fig. 4. The modes clearly rotate at 
different speeds and their rotation rates increase with mag- 
netic field strength. By comparing the rotation rates with those 
of higher azimuthal harmonics obtained from the high den- 
sity midplane coil array (seen as the series of dots at z = 
in Fig. 2) we find that the harmonics are not phase locked. 
Hence, the signal is not likely due to a passing vortex as was 
observed in Sisan et al. [10]. They are also not the result of 
the MRI since the least-damped mode should be axisymmet- 
ric. 

A least squares fit of the observed mode rotation rates to 
the real frequencies of the fast and slow MC waves gives an 
estimate of the local wave vector components and the fluid ro- 
tation rate for a given shear. The wavenumber components 
are fit since the observed magnetic field is much smoother 
than the velocity field due to the low Pm and cannot reveal 



detailed structures like the boundary layers. The fit is insensi- 
tive to the vorticity parameter £ and so we assume it to be zero 
for marginal stability consistent with hydrodynamic observa- 
tions [23]. The fit parameters result in a mode with k z = n/2h 
and k r = 71 / fa — r\) which corresponds to the smallest (and 
therefore least damped) radial wavenumber that can fit in the 
radial gap and a quarter-wavelength vertically. The growth 
rate determined from the fit is shown in Fig. 4b. Aside from 
an expected Doppler shift for the nonaxisymmetric modes, the 
MC wave model provides an excellent fit to the observations. 
From the growth rate, we find that the slow wave has a small 
positive growth rate for 1 kG but is otherwise damped. The 
model suggests that the MRI growth rate is too small at the 
rotation rates achieved to observe it with our diagnostics. By 
adjusting the model parameters we can predict the necessary 
rotation rate and magnetic field strength to observe the MRI 
based on the empirical observations of damped waves. 

Nonaxisymmetric waves have been observed in the 
PROMISE magnetized Taylor-Couette experiment [24] which 
has also observed the helical MRI [11, 25]. Sisan et al. [26] 
also observed rotating nonaxisymmetric spherical harmonic 
patterns in a turbulent liquid sodium spherical Couette flow 
which they attribute to the MRI, but may be due to a shear in- 
stability of the secondary meridional flow [27]. Runout of the 
inner cylinder could lead to a spinover mode due to the ellip- 
tical instability, but it is suppressed by the magnetic field [28]. 
It is still unclear why these nonaxisymmetric waves are fa- 
vored over the axisymmetric ones, and its resolution may re- 
quire 3D simulations as well as internal flow and magnetic 
field measurements. 

One conjecture regarding the source of these damped waves 
is that the turbulent flow provides perturbations of a broad 
range of wavelengths, but that the geometry of the vessel dic- 
tates which modes are realized [29]. Such is observed when 
a precessing top cap is used to drive inertial waves in a cylin- 
der filled with water [30]. The observed waves in this case 
are cavity resonances driven by unstable flow at the bound- 
ary discontinuities and not an instability of the bulk flow. It 
is also possible that these damped waves are driven through 
nonlinear wave coupling [31]. If such is the case, then these 
damped waves may be a saturation mechanism for the MRI. 
Although we have not observed the MRI in the experiment, 
nonlinear simulations of single mode MRI found that satura- 
tion was achieved by amplifying the vertical field through an 
a effect [32]. It is due to the ambiguity of the source of these 
waves that we are continuing to pursue an observation of mag- 
netically induced instability in a hydrodynamically quiescent 
flow. 

In summary, we have observed rotating modes in a turbu- 
lent Taylor-Couette flow of liquid metal that we identify as 
the fast and slow MC wave. We have identified a relationship 
between the slow MC wave and the MRI and have proposed 
a method of determining the threshold for instability through 
observation of driven MC waves. MC waves will be impor- 
tant in identifying the MRI in further experiments and may 
also play a role in saturation of MHD turbulence in rotating 
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FIG. 3: (a) Magnetic field measured by a single Mirnov coil, (b) The time series of (0,1) and (1,1) mode amplitudes from a 2D Fourier 
decomposition of the radial magnetic field where the notation (n,m) describes the vertical mode number n and the azimuthal mode number m. 
(c) The time series of mode phases. 
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FIG. 4: The real frequency (a) and growth rate (b) determined by 
fitting the phase speed of the two nonaxisymmetric modes observed 
for a range of applied magnetic field strengths to Eq. 1 . Error bars 
reflect the uncertainty in the linear fit to the phase as a function of 
time. The dashed red and solid blue lines show the fit of the fast 
and slow MC waves from Eq. 1 with values of £ = 0, k = 0.246 ± 
0.001 cm -1 , 6 = 1.336 ±0.007 rad where k z = £cosfl, and a = 
5.4 ± 0.9 rad/s. The shaded areas express the uncertainty in the fit. 
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